Anomalous Monopole In an Interacting Boson System 
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Anomalous monopole of disk shape is found to exist in the semiclassical theory of a two-mode inter- 
acting boson system. The quantum origin of this anomaly is the collapsing or bundling of field lines 
of Berry curvature caused by the interaction between bosons in the semiclassical limit. The signifi- 
cance of this anomalous monopole is twofold: (1) it signals the failure of the von Neumann- Wigner 
theorem in the semiclassical limit; (2) it indicates a breakdown of the correspondence principle 
between quantum and classical dynamics. 
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Magnetic monopole was first discussed by Dirac as a 
quantization condition for electric charge [J- Although 
it has never been observed in experiment as an funda- 
mental particle, the monopole has fascinated physicists 
ever since [2[. Interestingly, monopoles have attracted 
great attention in a very different context as degenera- 
cies or diabolical points of energy levels in parameter 
space@, Si- 

The examples include the degeneracy of 
Bloch bands in the Brillouin zone[6[ and energy levels in 
molecular magnets [7, 8]. The monopoles in this context 
are found to be crucial to understanding these systems. 

In this work we study a two-mode interacting boson 
system that depends on three external parameters. For 
simplicity, we focus on its ground state, which is doubly 
degenerate at one isolated point in the parameter space. 
We find that the field lines of Berry curvature emanating 
from the point are curved due to the interaction between 
bosons. Moreover, at large N limit, that is, when the 
number of bosons N increases to infinity, the field lines 
collapse and bundle into a two-dimensional disk whose 
radius is determined by the interaction strength. 

At large N limit, this boson system can be well de- 
scribed by a mean- field theory d, [Io|. Since this boson 
system belongs to a class of quantum systems which be- 
come classical at large N limit [ll|, this mean-field can 
be regarded as a semiclassical theory. We discover that 
the semiclassical (or mean-field) ground state of this bo- 
son system is degenerate at every point on the two- 
dimensional disk mentioned above. This means that the 
whole disk is a monopole. This is in stark contrast 
with what is demanded by the von Neumann- Wigner 
theorem (T3|: the monopole in a three-dimensional pa- 
rameter space is always a point-like object. Therefore, 
this anomalous monopole of disk shape indicates that 
the von Neumann- Wigner theorem fails in the semiclas- 
sical limit. Our further analysis shows that the magnetic 
charge is not uniformly distributed in the disk while its 
total charge is still 2ir, the Chern number In addition, 
this anomalous monopole is compared to an anomalous 



monopole that is formed in a trivial fashion. 

The Berry curvatures are computed for this system 
within the semiclassical theory and compared to the re- 
sults in the quantum description. The matching be- 
comes better as N increases as expected from the cor- 
respondence principle between quantum and classical 
dynamics [l3j . However, on the monopole disk, the Berry 
curvature differs significantly between its semiclassical re- 
sult and quantum result even in the large N limit. This 
shows that the existence of the anomalous monopole indi- 
cates a breakdown of the correspondence principle. This 
breakdown is analyzed from a fresh perspective by re- 
garding the three external parameters as the dynamical 
variables of a massive classical particle. 

The two-mode interacting boson system is described 
by the following second quantized Hamiltonian 



H N = y(a t 6 + aS t ) + ^(afe t - a r b) + 

A 



■(tfa-tfb) (tfa-tfb) 2 , 



AV 
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where a' , a and w , b are bosonic operators for two differ- 
ent quantum states, respectively, A > is the interacting 
strength between bosons, and V is the volume of the sys- 
tem. The three parameters, X,Y, and Z, characterize 
the influence from environment or another system. This 
Hamiltonian has its root in modelin g th e Bose-Einstein 
condensates in a double- well potential[irj . It also belongs 
to a class of Hamiltonians studied in Refs. @, H| for sin- 
gle molecule magnet if we introduce J x — {a)b + b'a)/2, 
J y = i(a>b - tia)/2^ and J z = (tfa - b%)/2. The focus 
of studies in Refs.0, Q is on the pattern and topologi- 
cal property of the monopoles. In this work we examine 
the "magnetic" fields, i.e., Berry curvatures, generated 
by the monopoles and their behavior in the semiclassical 
limit N — ► oo. Note that large N limit is always taken 
by keeping N/V constant. 

For simplicity, we concentrate on the ground state of 
this system. At point X = Y = Z = 0, we have Hn = 



2 




FIG. 1: Field lines of Berry curvature near the monopole 
for (a) the second-quantized Hamiltonian (dashed lines are 
for N = 2 and solid lines are for N = 10); (b) the mean- 
field Hamiltonian. Due to the symmetry around Z-axis, the 
Y component is omitted, c = 1. The lines in (b) are not 
straight as they appear. 



— A(a^a — wb) 2 /4V", whose ground state is either (a' a) = 
N or (b'b) = N. This means that the ground state of this 
boson system is doubly degenerate at point X = Y = 
Z = 0. The ground state is not degenerate elsewhere in 
the parameter space. 

This degenerate point or monopole &tX = Y = Z = 
generates "magnetic" field Bjv (Berry curvature) in the 
parameter space spanned by X, Y, Z. As one usually uses 
field lines to illustrate a magnetic field, we have computed 
numerically the field lines of Berry curvature and plotted 
them in FigQJa). For clarity, only the results for N = 2 
and N = 10 are plotted. Nevertheless an interesting 
trend is clearly demonstrated: the field lines are curved 
towards a disk defined by V X 2 + Y 2 = c = NX/V and 
Z = 0; the curving gets stronger as N increases. In 
fact, our numerical results show that the field lines will 
collapse and bundle (or converge) into the disk when TV 
approaches infinity. As we know, a magnetic monopole 



(or an electric charge) can be viewed as the converging 
point or the emitting source of field lines. This collapsing 
(or converging) of field lines suggests that the whole disk 
become a monopole in the limit N — > oo. This is indeed 
the case as we shall show. 

At large N limit, this boson system becomes "classi- 
cal" and can be described by the following mean-field (or 
semiclassical) Hamiltonian 



w tne i< 



H K = 



lim — — = — (a* b 

N^oo N 2 
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where a and b are complex amplitudes for the system in 
the two quantum modes. The normalization is one, i.e., 
| a, | 2 + \b\ 2 = 1. This kind of nonlinear Hamiltonian also 
appears in photoassociation systems [13. Il5j]. 

Within this semiclassical description, the ground state 
of this system is given by 




where p is the solution of the following equation, 



pVX 2 + Y 2 = {Z + cp)^l^ 



p 2 . 



(3) 



(4) 



This equation has one real root when v 1 X 2 + Y 2 > c. 
When X 2 + Y 2 < c 2 , it can have three real roots. In 
0, two of the three real roots given 
Y 2 )/c 2 have the same energy and 



particular, when Z 
byp= ±y/\ - {X 2 
are for the ground states. This means that in the semi- 
classical description of the system, the ground state is 
degenerate on the disk given by X 2 + Y 2 < c 2 and Z = 0. 
In other words, the whole disk is a monopole(see Fig(2]). 




FIG. 2: Anomalous monopoles. The disk is the monopole for 
the ground state and the thick vertical line is the monopole 
for the highest eigenstate. 

This anomalous disk-shaped monopole is very surpris- 
ing. According to the von Neumann- Wigner theorem, 



3 



the accidental degeneracy of a quantum system occurs 
only at isolated points in a three dimensional parameter 
space [l2|]. In other words, the monopole must be a point. 
This is indeed the case when our system is described by 
the second-quantized Hamiltonian: the ground state is 
degenerate only at point X = Y = Z = as we have 
already pointed out. However, in the semiclassical de- 
scription, that is, at large N limit, the monopole is a 
two-dimensional disk as the result of the curving and 
collapsing of field lines. This shows by example that the 
von Neumann- Wigner theorem does not hold in the semi- 
classical limit. 

The significance of this anomalous monopole can be 
appreciated from a different angle. We consider the high- 
est eigenstate of this system. For this eigenstate, its semi- 
classical monopole is a line as shown in Fig (2] However, 
this line-shaped monopole which appears also anomalous 
is trivial, not as significant as the disk-shaped monopole. 
The reason is as follows. In the second quantized model, 
the highest eigenstate has N equally spaced degenerate 
points along Z axis between — c and c. As N approaches 
infinity, these degenerate points merge into a line, form- 
ing in a trivial fashion the one dimensional monopole. 

Let us now examine this disk-shaped anomalous 
monopole in detail. Although the semiclassical Hamil- 
tonian H s is nonlinear, the Berry curvature 3 of this 
monopole can be computed as in a linear system [l6j|. 
That is to compute the curl of the vector potential 
A = (<^|V|<^>) with \(j>) given in Eq.(j3|. The Berry curva- 
ture 3 is found to be 



3 = 



P 



2{cp+ Z) 2 (cp 3 + Z 



■(R + cpz), 



(5) 



where R = {X, Y, Z} and z is the unit vector for Z direc- 
tion. This result is plotted as field lines in Figfljb). It is 
apparent that these semiclassical field lines away from the 
monopole disk are very similar to the field lines obtained 
with the second quantized model. Note that 3 has two 
different values on the monopole disk due to the double 
degeneracy of the ground state. By integrating 3 over 
a closed surface around a small area in the disk, we find 
that the "magnetic" charge is not uniformly distributed 
over the disk. The charge distribution is 



f 



cVc 2 -(A 2 + Y 2 ) 



(6) 



The integration of this charge density over the whole 
disk gives us a Chern number of 2tt. So, although the 
monopole has changed from a point to a disk as the 
semiclassical limit is approached, the total charge does 
not change. It is worthwhile to mention that the total 
charge of the line-shaped monopole is infinite as easily 
inferred from its trivial origin. 

Berry [l3} once established a semiclassical relation be- 
tween Berry phaseQ and Hannay's angle 17, [l8[ in ac- 



cordance with the correspondence principle. This semi- 
classical relation basically says that the two-forms for 
Berry phase and Hannay's angle (the two-form for Berry 
phase is the usual Berry curvature) are the same in 
the semiclassical limit % ~ * 0. This semiclassical rela- 
tion should hold in this interacting boson system. We 
define two pairs of conjugate variables, p a — s/iha*, 
q a = Viha and pb = Vihb* , qb = Vihb for the semi- 
classical Hamiltonian [l9l| . The quantization is realized 
with the following commutators, 



[q a ,Pa] = [qb,Pb] = ih/N. 



(J) 



One can obtain the second quantized Hamiltonian([T]) 
with the following substitution a = y/ N/ihq a , a* = 
y/N/ihp a and b = y/N/ihq b , = y/N/ihp b . These 
commutators show why N — > oo is the semiclassical 
limit. As a result, the semiclassical relation established 
by Berrydl for this boson system is 



lim SB = lim (— - 3) = 0, 



(8) 



Note that the Hannay's angle in the semiclassical sys- 
tem H s is just the Berry phase generalized for nonlinear 
quantum system in Ref. [161] and the semiclassical Berry 
curvature 3 is the two-form for this Hannay's angle. 

Our numerical results show that the relation ([5]) indeed 
holds almost everywhere in the parameter space except 
on the monopole disk. On the disk, the semiclassical 
Berry curvature 3 has a non-zero I component while the 
quantum Bjv always points radially in the Z — plane. 
Furthermore, the quantum Berry curvature Bjv diverges 
exponentially with N on the monopole disk while the in- 
plane component of the semiclassical 3 does not. We 
define d = \SB |, where the superscript I denotes the 
component of the vector parallel to the XY plane. The 
difference d is plotted in Fig [31 where we see the differ- 
ence d increases exponentially with N. This diverging 
difference shows that the semiclassical relation in Eq.© 
is broken. Therefore, the disk-shaped monopole also sig- 
nifies the breakdown of the corresponding principle be- 
tween quantum mechanics and classical mechanics. In 
the following, we shall look into this breakdown from a 
very different angle, showing that some quantum effect 
remains in the semiclassical limit. 

We treat the three parameters, R = {X,Y,Z}, as 
dynamical variables of a massive and classical particle, 
whose Hamiltonian is H c — P 2 /2M+ t/(R). In this way 
we obtain a Born-Oppenheimer type system where a fast 
quantum system is coupled to a heavy and slow classical 
system [19( 



(9) 



where \^) is a general wavevector for the boson system. 
It is reasonable to expect that at the large N limit we 
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FIG. 3: The difference d = \8~B l \ as a function of the number 
of bosons N at a point on the monopole disk. The inset shows 
the results for points away from the disk, c = 1. 



can simplify the above Hamiltonian by replacing its first 
part by the semiclassical Hamiltonian in Eq.([2]), 



H' = NH S + H c 



(10) 



However, in the following, we shall show that the two 
Hamiltonians H and H' do not have the same dynamics 
even in the limit TV — > oo. That is, there is a difference 
between the two Hamiltonians which does not vanish as 
N increases. 

For a Born-Oppenheimer system, it is well known that 
the dynamics of the slow system H c is controlled by two 
forces. Born-Oppenheimer force and geometric forceps. 
M,^M- The Born-Oppenheimer force F BO = -V£„(R) 
is the gradient of an eigenenergy of the fast system. If B 5 
is the Berry curvature of the fast system, the dynamics 
of the slow system is then given by 



MR 



' BO 



RxB 9 - VC/(R) . 



(11) 



One can use either the second quantized Hamiltonian 
i/jv or the semiclassical Hamiltonian H s to compute both 
forces. 

We consider a special case, where the slow classical par- 
ticle is set with the initial condition, X = r < c, Y = 0, 
Z = 0, and X = 0, Y = v and Z = while the fast boson 
system is kept in its ground state. We set up ourselves to 
a task to design a potential U(R) so that the slow par- 
ticle stays in Z — plane and makes a circular motion. 
When N is large, one would feel safe to design f7(R) by 
using the semiclassical Hamiltonian H' to compute F bo 
and B 5 . However, due to the exponential breakdown of 
Eq.© discussed above, such designed U(R) will not be 
able to keep the classical particle in the Z = plane: 
the strong parallel component of Bat will kick the par- 
ticle out of the plane. This shows that there is always 
some physical consequence which can not be counted in 
the semiclassical theory. It is interesting note that Berry 



and Robbins once pointed out that a kind of friction in a 
chaotic classical system does not exist in its correspond- 
ing quantum system and called it discordance [2lj]. What 
we observe here is similar to this discordance although 
our system is not chaotic in the semiclassical limit. 

In conclusion, we have found an anomalous monopole 
of disk shape in a two-mode interacting boson system. 
This kind of anomalous monopole should exit in a gen- 
eral interacting boson system. For example, if one man- 
ages to compute the Bloch bands of an interacting bo- 
son system in a three dimensional periodic potential, 
one should expect such an anomalous monopole in the 
Brillouin zone. We have further demonstrated that this 
anomalous monopole is an indication of the failure of the 
von Neumann- Wigner theorem in the semiclassical limit 
and the breakdown of the correspondence principle be- 
tween quantum and classical dynamics. 
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